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Abstract 

We study the role of the noise in the dynamics of two competing species. We con- 
sider generahzed Lotka-Volterra equations in the presence of a multiphcative noise, 
which models the interaction between the species and the environment. The inter- 
action parameter between the species is a random process which obeys a stochastic 
differential equation with a generalized bistable potential in the presence of a peri- 
odic driving term, which accounts for the environment temperature variation. We 
find noise-induced periodic oscillations of the species concentrations and stochastic 
resonance phenomenon. We find also a nonmonotonic behavior of the mean extinc- 
tion time of one of the two competing species as a function of the additive noise 
intensity. 

Key words: Statistical Mechanics, Population Dynamics, Noise-induced effects 
PACS: 05.40-a,87.23Cc,89.75-k 



1 Introduction 



Recently the influence of noise on population dynamics has been the subject 
of intense theoretical investigations [1,2,3,4,5]. The problem of the stability 
of the biological complex systems in the presence of noise for example is one 
of the more discussed in research activity of these last years [6]. The noise 
through its interaction with the nonlinearity of the systems can give rise to 
new counterintuitive phenomena like stochastic resonance [7], noise enhanced 
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stability [8], noise delayed extinction, etc... The principal aim of this work 
is the investigation of a generalized Lotka-Volterra model with a random in- 
teraction parameter between two competing species and in the presence of a 
multiplicative noise [1]. We find that the noise has a constructive role. In fact 
it is responsible for the generation of quasi-periodic temporal oscillations and 
the enhancement of the response of the system to a driving force producing 
stochastic resonance. Moreover we find noise delayed extinction, i. e. a non- 
monotonic behavior of the average extinction time of one of the two species 
as a function of the noise intensity. 



2 The model 



Time evolution of two competing species is obtained within the formalism of 
the Lotka-Volterra equations [9] in the presence of a multiplicative noise 



dx 

— = fiix{ai-x- f3i{t)y) + x^^{t) 
dy 

— = H2y{a2-y - f32{t)x) + yiy{t), 



(1) 

(2) 



where ^xif) and ^y{t) are statistically independent Gaussian white noises with 
zero mean and correlation function {Ci(t)^j(t')) = a6(t — t')6ij = x, y). The 
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Fig. 1. The bistable potential U {(3) of the interaction parameter (3 {t) . The potential 
U{(3) is centered on (5 = 0.99. The parameters of the potential are h = 6.25 • 10~^, 
= 0.05, p = -0.01. 



time series for the two populations are obtained setting ai = a2 = a, = 
P2it) = f^if)- It is known that for /3 < 1 a coexistence regime takes place, that 
is both species survives, while for /5 > 1 an exclusion regime is established. 
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that is one of the two species vanishes after a certain time. Coexistence and 
exclusion of one of the two species correspond to stable states of the Lotka- 
Volterra's deterministic model [10]. Real ecosystems are immersed in a noisy 
nonstationary environment, so also the interaction parameter is affected by the 
noise and some other deterministic periodical driving such as the temperature. 
The change in the competition rate between exclusion and coexistence e. g. 
occurs randomly because of the coupling between the limiting resources and 
the noisy environment. A random variation of limiting resources produces a 
random competition between the species. The noise therefore together with 
the periodic force determines the crossing from a dynamical regime (/3 < 1, 
coexistence) to the other one {P > 1, exclusion). To describe this continuous 
and noisy behaviour of the interaction parameter (3{t) we consider a stochastic 
differential equation with a bistable potential and a periodical driving force 

+ 5,(0. (3) 



where U{P) is a bistable potential (see Fig.l) 

Ui/3) = hi/3 - (1 + p)r/v' - 2hif3 - (1 + p))Vr/^ (4) 

and h is the height of the potential barrier. The periodic term takes into ac- 
count for the environment temperature variation. Here 7 = 10~^ and co'o/(2vr) = 
10~^. In Eq.(3) ^p{t) is a Gaussian white noise with the usual statistical prop- 
erties: i^fsit)) = and (^/3(t)^/3(t')) = (^H^ - ^')- Due to the shape of U{(3) it 
is reasonable to expect a coexistence regime for j3{0) < 1, when deterministic 
case {(.(sit) = 0) is considered. 

2.1 Stochastic resonance 

First we investigate the effect of the noise on the time behavior of the species. 
Since the dynamics of the species strongly depends on the value of the inter- 
action parameter, we initially analyze the time evolution of (3{t) for different 
levels of the additive noise a p. Specifically for ct/j = we obtain a periodical 
behavior of l3{t) in the coexistence region (see Fig. 2a). For low noise intensity 
{(Tp <^ h) we get the same periodical behavior of deterministic case, slightly 
perturbed by the noise (Fig.2b). For higher noise intensity (a^ ~ h) the pe- 
riodical behavior of the interaction parameter jumps between the two values 
/3 = 0.94 < 1 and (3 = 1.04 > 1, which characterize the coexistence and the 
exclusion regimes. We see in Fig. (2c) the typical picture of stochastic reso- 
nance. A further increase of the noise intensity produces a loss of coherence 
and the dynamical behavior is strongly controlled by the noise (Fig.2d). There 
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Fig. 2. Time evolution of the interaction parameter for different values of the 
additive noise ap. (a) up = 0; (h) = 1.78 • 10"^; (c) = 1.78 • 10"^; (d) 
= 1.78 • 10^^. The values of the parameters are: 7 = 10""*^, ujq /{2t:) = 10""^. 



is therefore a cooperation between periodical driving of the temperature, due 
to some geological cause, and the environmental noise. The noise can synchro- 
nize with the periodical driving giving rise to the stochastic resonance effect, 
which affect considerably the dynamics of the ecosystem. The exact value of 
= 1.78 ■ 10^'^ used in Fig. 2c is obtained from the statistical synchronization 
typical of the SR phenomenon 

Tk = To/2 (5) 



where Tk is the Kramers time given by 
27r 

Tk = I exp [2h/ ap] . (6) 

|[/"(0.99)|[/"(0.94) 



and To is the period of the driving force. In Eq.(6) f/"(0.99) and f/"(0.94) 
are the second derivative respectively calculated in the unstable and stable 
states of the potential. To analyze the dynamics of the two species we fix 
the additive noise intensity at the value cr^ = 1.78 ■ 10~^, corresponding to 
a competition regime between the two species periodically switched from co- 
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Fig. 3. Time evolution of both populations at different levels of the multiplicative 
noise: (a) a = 0; (b) a = 10-^^• (c) a = 10-i°; (d) a = IQ-^; (e) a = lO""^; (f) 
a = 10^^. The values of the parameters are /x = 1, a = 1, 7 = lO""*^, ujq/2tt = 10~^. 
The intensity of the additive noise is fixed at the value ap = 1.78 • 10~^. The initial 
values of the two species are 2;(0) = y(0) = 1. 



existence to exclusion. The temporal series of the two species are obtained 
for different values of the multiplicative noise intensity a = = <Jy The 
initial values of the two species are x(0) = y(0) = 1. After a short transient 
in the coexistence regime, which is not visible in Fig. 3 because of the scale 
used for the x axes, both species reach the value: Xst = = + ~ 1/2- 
From then on the densities of both species oscillate quasi regularly around 
this value with amplitudes dependent on the multiplicative noise intensity. In 
particular for a = (see Fig. 3a) and for very low levels of multiplicative noise 
(a = 10~^^) the two species coexist and the populations show correlated oscil- 
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lations. For higher level of the multiplicative noise {a = 10"'*) the amplitude 
of the oscillations increases, anti-correlated oscillations appear and random 
periodical inversions of populations occur (see Fig. 3c). A further increase of 
the noise {a = 10~^) produces a degradation of the signal and a loss of co- 
herence of the temporal series for the species (see Fig. 3d). In the absence of 
multiplicative noise the Lotka-Volterra equations are symmetric. The choice 
of the same values for the initial conditions of the two species densities main- 
tains this symmetry, so the species oscillate in phase around the stationary 
value a/{l + (3), even if the additive noise that controls the SR effect of the 
interaction parameter was fixed at a value that produces many transitions to 
the exclusion regime (see Fig.3a). In the presence of very weak multiplicative 
noise (cr = 10~^^) the species oscillate in phase almost regularly with some 
little anticorrelated spikes. By increasing the amplitude of the multiplicative 
noise we obtain an enhancement of the amplitude of both species, as we can 
see from Figs. 3b, 3c and 3d. These figures show this transient behaviour from 
a regime of quasi correlated oscillations (Fig. 3b) to a regime of anticorrelated 
oscillations (Fig.3d). The presence of the multiplicative noise breaks the sym- 
metric dynamical behaviour of the ecosystem. For /5 > 1, i. e. in the exclusion 
regime, this symmetry breaking determines the "divergent" behaviour of the 
trajectories of the two species: one tends to survive, while the other one tends 
to extinguish. From the inspection of Fig. 3 we note that in the presence of 
a bistable potential modulated by a weak periodic force the response of the 
system may be enhanced by the presence of the noise. The periodicity of 
the noise-induced oscillations in the time behavior of the population densities 
shown in Fig.3e is the same of the driving periodic term of Eq.(3). This is 
the signature of stochastic resonance (SR). In order to underline the presence 
of SR we analyze the squared difference of population densities {x — y^. In 
Fig. 4 it is shown the SNR of this quantity as a function of the multiplicative 
noise intensity a, for = 1.78 ■ 10~^. We calculate the SNR of {x — yY be- 
cause of the random periodical inversion of population densities. We note that 
dynamics of {x — y) is mainly affected by the multiplicative noise, as we can 
see from Eqs.(2). The SNR has been obtained by performing 190 realizations 
of Eqs. (2). A maximum at a = 10~^ is present. From the above analysis it 
is clear the role of the two noise sources: the additive noise determines the 
conditions for the different dynamical regimes of the two species, the multi- 
plicative noise produces a coherent response of the system by a mechanism of 
symmetry breaking of the dynamical evolution of the ecosystem. 

2.2 Noise delayed extinction 

Now we investigate the average extinction time of one species as a function of 
the noise intensity ci/j. To do this we fix the multiplicative noise at a low level 
intensity in such a way that the system is far enough from the SR regime and 
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Fig. 4. Log-Log plot of SNR as a function of the multiplicative noise intensity. 
The SNR corresponds to the squared difference of population densities (x — y)'^. The 
values of the parameters are the same of Fig. (3). 

the dynamics is not strongly perturbed by the noise. So we choose a = 10~^ 
because we are not interested in the coherent behavior of the ecological system 
but we are focused on the effect of the additive noise on the average extinction 
time of the species. We do simulations of Eqs. (1) and (2) by performing 200 
realizations. The results are shown in Fig. (5). We use the bistable potential 
U{l3) of Fig.l with the initial condition /3(0) = 0.94. In this condition the 
ecosystem is in the coexistence regime, that is the deterministic extinction 
time of both species is infinite. By introducing noise, exclusion takes place 
and a finite mean extinction time (MET) appears. By varying the intensity of 
the additive noise in Eq.3 we obtain, of course, a variation of the average ex- 
tinction time. The delayed extinction is obtained for noise intensities ranging 
from the intermediate regime (2 in Fig. 5a) to the coexistence regime obtained 
with higher values of a (3 in Fig. 5a). In real ecosystems we have always a 
given noise intensity, which corresponds to a finite mean extinction time. Due 
to some environmental cause the noise intensity can considerably change, as 
it is observed in experimental data of populations in a very long time inter- 
val [12]. Therefore the dynamical behaviour shown in Fig. 5 should explain such 
physical situations, where the variation of the environmental noise produces 
a delayed extinction of some population. By increasing the noise intensity we 
obtain noise delayed extinction and the average extinction time grows reaching 
a saturation value, which corresponds to a situation where the potential bar- 
rier is absent. We find nonmonotonic behaviour of the MET as a function of 
the noise intensity a^, with a minimum value Tmin = 40.47 at = 2.75 ■ 10~^, 
which is of the same order of magnitude of the barrier height h (see Fig. 5a). 
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Fig. 5. ('aj Mean extinction time of one species as a function of the noise inten- 
sity ap. Time evolution of both species for different levels of additive noise: (h) 
Gjs = 10~^, (c) = 2 • 10"'^, (d) Gjs = . The values of the parameters are 
fj, = 1, a = 45, 7 = 10^^, a;o/27r = 10"'^. The intensity of the multiplicative noise is 
fixed at the value a = 10~^. The initial values of the two species are x{0) = y{0) = 1. 

The Kramers time corresponding to this noise intensity is = 41.6, that is 
approximately equal to Tmm- This result is due to the noise driven dynamics. 
In fact for low value of the noise intensities the average time to overcome the 
potential barrier is very high, i.e. long Kramers times. The ecosystem remains 
in the coexistence regime for a long time and the extinction time is very large. 
For noise intensity of the same order of magnitude of the barrier height the 
system goes towards the exclusion regime of one of two species and the average 
extinction time is approximately equal to the Kramers time. We get the min- 
imum value of MET. For higher values of noise intensity the Kramers time 
becomes very small and the representative point of the (3 parameter moves 
between the two minima in a very short time. In this condition the system 
"sees" the average value of the interaction parameter (/? = 0.99), which gives a 
coexistence regime. In Figs.5b,c,d we show the time evolution of the ecosystem 
corresponding to the points 1, 2 and 3 of Fig.Sa. We have coexistence regime 
in points 1 and 3 and an exclusion regime in point 2. We report finally in the 
following Fig. 6 the behavior of the average extinction time as a function of the 
Kramers time. This figure confirms our physical description of the extinction 
dynamics of the system. 
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Fig. 6. Semilog plot of the mean extinction time of one species as a function of the 
Kramers time. 

3 Conclusions 



We report a study on the role of the noise in the dynamics of two competing 
species. We consider two noise sources: a multiplicative noise and an additive 
noise, which produces a random interaction parameter between the species. 
The noise induces a coherent time behavior of two species giving rise to tem- 
poral oscillations and enhancement of the response of the system through 
stochastic resonance phenomenon. Specifically the additive noise controls the 
switching between the coexistence and the exclusion dynamical regimes, the 
multiplicative noise is responsible for coherent oscillations of the two species. 
The SR in the dynamics of interaction parameter (3 induces SR phenomenon 
in two competing species. These time behaviors are absent in the deterministic 
dynamics. We note that our model is useful to describe physical situations in 
which the amplitude of the periodical driving force, due to the temperature 
variations, is weak and therefore unable to produce considerable variations of 
the dynamical regime of the ecosystem. The synergetic cooperation between 
the nonlinearity of the system and the random and periodical environmental 
driving forces produces therefore a coherent time behaviour of the ecosystem 
investigated. The noise is also responsible for a delayed extinction which gives 
rise to a nonmonotonic behavior of the average extinction time as a function 
of the additive noise intensity. We note finally that these noise induced effects 
should be useful to explain the time evolution of ecological species, whose 
dynamics is strongly affected by the noisy environment [5,6,11,12]. 
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